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Abstract 

We use a macroscopic Hamiltonian approach to study the pinning of a solid-liquid-vapour 
contact line on an array of equidistant stripes of obstacles perpendicular to the liquid. We 
propose an estimate of the density of pinning stripes for which collective pinning of the 
contact line happens. This estimate is shown to be in good agreement with Langevin equation 
simulation of the macroscopic Hamiltonian. Finally we introduce a 2-dimensional mean field 
theory which for small strength of the pinning stripes and for small capillary length gives an 
excellent description of the averaged height of the contact line. 
PACS numbers: 68.10.-m, 68.45.Gd,68.45Ws. 



1 Introduction 



The spreading (wetting) of a liquid on a solid is important in a widespread field of prac- 
tical applications such as lubrication, the efficiency of detergents, oil recovery in a porous 
medium and the stability of paint coatings 0. The motion of the interface is often extremely 
sensitive to impurities and roughness, which tend to pin (stick) the interface. Different situ- 
ations arises for the motion of the boundary line between a solid, liquid and vapour (called 
the triple line or the contact line) depending on the heterogeneousity of the solid and de- 
pending on whether the liquid completely wets the solid surface or incompletely wets the 
solid [3- 13]. Due to the presence of a microscopic precursor film that advances ahead of the 
macroscopic liquid, the case of complete wetting is noticeable insensitive to the heteroge- 
neousity of the solid. In the incomplete wetting case obstacles tend to pin the contact line 
which makes the statics and dynamics of the contact line highly sensitive to the specific 
form of the heterogeneousity of the solid. Similar kinds of problems are also encountered 
in other situations where an elastic body is pinned by a random potential |14[ such as flux 



pinning of type-II high T c superconductors [HJ, pinning of charge density waves p6| , pinning 
of magnetic domain walls [IT] and dislocation pinning [18]. 



We have proposed a macroscopic Hamiltonian approach to study the pinning of a solid- 
liquid-vapour interface when a solid is pulled vertically out of a liquid (called the immersion 
geometry). We investigate the simple case of the statics of an array of equidistant stripes of 
obstacles, and we study how the averaged height of the triple line, depends on the capillary 
length and the density and strength of the stripes. We hope by studying this simple case 
to understand some of the important physics that is involved in the more complex case of 
completely random pinning sites. 

2 A macroscopic Hamiltonian approach to the pinning 
of solid— liquid— vapour interfaces 

The geometry of the problem is chosen so that the liquid-vapour interface is in the x-y- 
direction and the solid is pulled vertically out of the liquid in the ^-direction. The origin of 
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the z-axis is taken to coincide with the liquid-vapour interface. The local surface energies for 
the liquid-gas interface, ■jlg, the solid-liquid interface, jSL(y, z) and the solid-gas interface 
1sg{Vi z ) determine a local force balance expressed by Young's [|19| relation: 



iLsiViZ) +1LG COS e(y,z) = -isG{y,z) (1) 

6(y, z) is the local macroscopic contact angle between the solid and the liquid, 
The Hamiltonian reads 
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where p is the mass density of the liquid (assumed constant) and g the gravity constant. The 
first term in this expression is the gravity potential energy, the second term describes the 
surface energy between the liquid and the gas, and the last two terms account for the surface 



energies due to contact between the solid wall and the liquid, respectively, gas phase . 

If one considers how the Hamiltonian changes by a small change in 5h(x,y), the three 
parameters r )sL{.y^z) ) ^ SG {y ) z) and that enters the Hamiltonian, can be expressed in 
terms of just two variables, namely the liquid-gas surface energy, Jlg, and the macroscopic 
contact angle 9(y, h(x,y)) using Young's relation Eq. (jl]): 



pgh( Xl y) - 7lg V'/i(x, y)(l + (Vh(x, y)) 2 )^ 



5h(x,y) 

+ lLG V 2 h{x, y){l + (Vh(x, y)f)-'i{Vh{x, y)) 2 - lLG cos(%, h{x, y)))6(x\3) 

The sole quantity that has the dimensions of a length is called the capillary length and is 
defined as: 



On the length scales for which the Hamiltonian Eq. (0) is supposed to be valid, one can safely 
ignore fluctuations due to the temperature, so the problem of finding the configuration h(x, y) 
that minimizes Ti is a zero temperature problem. 
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One way to obtain the equilibrium configuration h(x, y) described by the Hamiltonian 
Eq. (HD is to perform simulated annealing using Monte Carlo simulation. Notice that since 
only the change of the total energy is needed in a Monte Carlo update, the equilibrium 
state is completely specified by the three variables pg, 9(y, z) and ^lg- In Ref. [[§] simulated 
annealing was carried out for the case 9(y, z) = constant, and for the case of equidistant 
stripes of obstacles (with contact angle 9') in the ^-direction. The case 9(y, z) = constant 
served as a check of the validity of Eq. (|2]) since the profiles h(x, y) can be directly compared 
to various analytical results^!]]. Initial findings for the case of stripes of obstacles indicated 
that the averaged height of the triple line, < h(0,y) >, was linear in the density of pinning 
sites, c, for small values of c and with a crossover to nonlinear behavior for c — > 1. The 
density for which the crossover happened, c*, was an increasing function of -. This is to 
be expected, since for small c a given pinning stripe does not feel the presence of the other 
pinning stripes, and < h(0,y) > can be obtained as a simple superposition over all the 
pinning stripes of the profiles of individual stripes that pin the liquid-vapour interface. On 
the other hand when c becomes larger the pinning stripes mutually (collectively) lift the 
liquid, and the resulting profile h(0,y) can not be obtained as a simple superposition over 
the pinning stripes. Therefore one expects that collective pinning sets in once the averaged 
distance between pinning stripes, d, becomes of the order of the capillary length k, giving 
c* ~ 1/d ~ 1/k. We will in this paper show that another length scale enters the problem, 
so that the before mentioned argument has to be modified. 

It turns out that simulated annealing has the disadvantage of requirering large amounts 
of computing time, since one needs to perform a slow annealing sequence r(t)p2|| in order 
to bypass metastable configurations. Furthermore the optimal sequence T(t) depends on 
the parameters 9, 9', c, k, the lattice constant a, and one should in principle determine the 
optimal sequence T{t) and the optimal step size Ah for each set of parameters values used 
in a simulation. In this paper we have instead used Langevin equation simulations of Eq. (|3|) 
since it was found to be computational more efficient than the Monte Carlo simulations. 
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That is, we have numerically solved the equation 

5t " 6h(x,y)' 

T is a mobility constant, t is the time which we let go to infinity in order to find the 
equilibrium state, and T is the free energy of the system which we will assume T ~ Ti,, 
since we neglect fluctuations due to the temperature. In order to take proper account of the 
boundary condition for the contact of the liquid with the solid plate one has that 

^^|_o=cot(%,MO,y))). (6) 

This condition has to be introduced via a Langrange multiplyer as an extra term to the 
Hamiltonian Eq. (|2]), or equivalently one leaves out the contact terms in the energy Eq. (|2|) 
(which amounts to leave out the term — ■jlg cos(9(y, h(x, y)))8(x) in Eq. ([$])) and expresses 
them instead via Eq. (|6])). 

When discretizing Eq. ([5]) the functional form of theta for a realization of stripes of 
obstacles takes the form 

9{V) z ) — ^modulus(y,i),O0' + (1 — ^modulus(j/,i),o)^) (7) 

with 9' the value for the contact angle on the stripes of obstacles, 9 the value of the contact 
angle between the stripes of obstacles, and 1 < i < L/a an integer determining the density 
(0 < c < 1) of pinning stripes. One notice that due to translational invariance in the z- 
direction, there is no pinning force acting when the liquid is pulled vertically out of the 
liquid, whereas a rotation of the solid about an axis in the x-?/-plane or a translation of the 
solid in the ^-direction produces a pinning force. 

3 Results 

In Fig. 1 is shown the averaged height of the triple line < h > versus the density c for different 
values of k from a numerical calculation of Eq. (|5],|5],|7]) . The discretization parameter was 
chosen a = 0.04 (the same for all the results represented in this paper), a typical time step 
At ~ 10 -4 — 10~ 2 depending on the values of k, 9, 9', c; the lattice sizes were chosen between 



Figure 1: The averaged height of the triple line < h > versus the density c for different 
values of k. k — 4(D), k = 2(+) and k = l(O) respectively. 6 = |, 6 1 ' = 1.5. 



Figure 2: The averaged height of the triple line < h > versus the density c for different values 
of k. k= l(O), « = 0.5(+), « = 0.25(D) and « = 0.125(x) respectively. = §,0' = 0.5. 
From Eq. ([□]) one finds: c curvatU re(K = 1-0) = 0.15, c curvatU rc(K = 0.5) = 0.30, c curvatU re(/« = 
0.25) = 0.60, and c curvaturc (fi: = 0.125) > 1. 



L/a x M/a = 100 x 100 and L/a x M/a = 100 x 400 with the larger lattice size for larger 
k and the number of time steps to reach equilibrium were in the range 2 x 10 4 — 10 6 . The 
distance between two pinning stripes, d, is given by the total length of the solid in the 
^-direction, divided by the total number of stripes: 

L a 

d=— = -. (8) 

cL/a c 

Therefore the critical density, c* , for which collective pinning should set in given by equating 
this distance d with the capillary length k is: 

< - - K P) 

For the given variables this gives the critical densities c* = 0.01, 0.02, 0.04. Since < h > 
increases linearly with c for all the three values of k without showing any sign of crossover, 
we conclude that we are in a single pinning regime contrary to the simple argument that we 
suggested above. 

The reason for this discrepancy is that the pinning stripes introduce a new length scale in 
the problem, namely the distance between the height to which the liquid would rise without 
pinning stripes < h > minus the height to which the liquid would rise with a density of one 
of pinning stripes, < h' >: 



H=<h' >~<h Q > = K(J(l-sm{0)-J{l-sm{0')). (10) 
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Therefore if H plays a role for the onset of the single/collective pinning regime one should 
be able to go from one pinning regime to another by changing either \9' — 9\ or k. This 
statement we have confirmed by keeping k constant and increasing the quantity \9' — 9\ for 
various values of k (see the discussion after Fig. 3). Furthermore, in Fig. 2 is shown the 
averaged height of the triple line < h > versus the density c for different values of k for a 
fixed value of \9' — 9\. One notice a crossover from a single pinning regime to a collective 
pinning regime as one increases c, and with the crossover appearing for smaller c the larger 
the value of k. For given values of k, 9, 9' and for small c's a stripe do not feel the presence 
of its neighbor stripes. As c increases, the curvature of the triple line between two stripes 
increases up to a point where the curvature becomes so large, that the cost in gravitational 
energy by lifting the liquid between two stripes is outbalanced by a decrease in curvature and 
thereby surface energy. When this happens the system is in the collective pinning regime. 
Assuming that the decay of the triple line away from a stripe is exponential, we estimate 
that curvature effects become important when the distance the triple line has decayed after 
one correlation length, becomes on the order of the averaged distance between two stripes: 

— — = d* = Hexp(-l) 

'-'curvature 

* a /-, -i ^ 

'-'curvature 1 1 ' ' 



- sin(0) - ^(1 - sm(9'))exp(-l) 
Assuming the form of the triple line in between two pinning stripes can be described as a 
segment of a circle, another way of stating Eq. flTTl ) is to say that collective pinning happens 
once the radius of curvature between the stripes becomes on the order of the averaged dis- 
tance between the stripes. In Fig. 2 and Fig. 3 are indicated the critical densities c* urvature 
for the various parameters used in the simulation. Giving that Eq. (|ll]) is only based on a 
simple order of magnitude argument, the agreement with the onset of the collective pinning 
deduced from the simulations is striking. Besides, using Eq. QTTj) for the simulations per- 
formed in Fig. 1, all give a c* urvature larger than one, meaning that collective pinning should 
not occur for these values of parameters, which is in agreement with Fig. 1. Based on these 



results, we conjecture that Eq. (|Tl| ) should also give a good estimate in the case of random 



pinning sites, with 9 describing the averaged value of the strengths of pinning sites, and 9' 



the fluctuations about it. 



4 Mean field theory 

We now propose a simple mean field theory in order to obtain the averaged height of the 
triple line, < ho(c, k, ao, 9, 9') >, as a function of the density of pinning stripes c. The mean 
field assumption amounts to consider just one pinning stripe, and to find the minimum 
energy of this configuration. Furthermore we simplify the problem, and consider only a 2- 
dimensional projection of the 3-dimensional liquid meniscus onto the solid; We will assume 
the functional form for the height h as a function of distance y away from the pinning stripe: 

y G [0 : a] : h(y) = h ; y G [a : d] : %) = h exp(-^) + h x (12) 

The two constants ho and /ii will be determined form the requirement that the mean field 
energy E mcmi of the system attains it minimum. E mcan is given by: 

]_ ra+d y 

-Emcan = T^Pd 1 ^ ~ Ilg cos(9')h a - / dy^ LG cos(9)h exp( — ) - >y LG cos(Q)hid 

Z J a K 

ra+d I 2y ra+d \ ra+d y 

+ dy-pgth exp( )+/ dy-pgth 1 + dypgthohexp^ ).(13) 

Ja Z K Ja Z Ja K 

The variable t is supposed to take into account the contribution in energy of the three 
dimensional meniscus on the two dimensional projection due to interfacial and gravitational 
energy. The form of t will be determined by the requirement: 

< h (c, K,a,9,9') h ,c — >■ ; < h (c, k, a, 9, 9') h 1: c — > 1 (14) 

Minimizing the energy 

dE dE 

^-^mean q_ WJ - / mean q 

dh ' dhi 

pgth\a - 7 LG cos(9')a(--f LG cos(6») + pgthi)K[exp(--) - exp(- a + - )] 

K K 

1 , . . 2a, . 2(a + d). 1 

+ -pgth K[exp( ) - exp( ) = 0; 

2 K K 

~1lg cos(9)d + pgtho) K[exp(— — ) — exp( -— )] + pgthid = 0. (15) 

K K 
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Eq. fll5| ) is two equations with two unknowns from which ho, hi can be determined. Now 
< ho{c, k, ao, 0, 9') > is given by the integral over the functional form Eq. fll2|) : 

pa pa+d y> 

(a + d) < ho(c, k, ao, 0,6') > = / dxh + / dx[hnexp( )+hi], (16) 

J Ja K 

which from the solutions of Eq. (|T5|) can be written: 

^ 2 1lg cos(0')ac 



< h Q (c, k, a , 0, 0') > 



- (ex P (-f ) - exp(-£))* + |«(exp(-f ) - exp(-%)] 

| 7 LG cos(g)(l-c) 



In order to fulfill the condition Eq. ([l^) it can be seen from Eq. ( |17D that a proper choice of 
t is 



* = \h—m ■ 77^-^ r-77^r(ccos(0') + (1 -c)cos(^)) (18) 

Y2p#[c(l-sin(fl') + (i-c)(l-sin(0)] v v ; v ; v is \ i 

Thus we end up with the form for < h > given by: 

. / fl fl /w k 2 7lg cos(#>c 

< finlC, k, a Q , 0,0 ) > = 5 — 

2^-^(exp(-^)- e xp(-^ + i«(ex P (-f)-e^(-^)] 
| 7 Lg cos(g)(l-c) 

Finally we will assume that the width of the pinning stripes a depends on the density of 
pinning stripes c in the following way: 

c E [0 : 0.5] :a = a ; c G [0.5 : 1] : a = a L2(c - 0.5) (20) 

which just states that the no two pinning stripes will be neighbors until c = 0.5 where after 
the width of a stripe increases linearally with c. 

In Fig. 3 is shown the averaged height of the triple line < h > versus the density c for 
different values of with a fixed value of k and 0' . One notice that the density for which 
collective pinning begins clearly depends on and is well described by Eq. ([TT|) . Furthermore 
the 2-dimensional mean field solution gives a good description of the 3-dimensional Langevin 
solution, with the best agreement for small H in Eq. flT0|). This is to be expected since for 
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Figure 3: The averaged height of the triple line < h > versus the density c for different 
values of 9. k = 0.5,6' = 1.0. 9 = 0.5(x), 9 = 1.2(A), 9 = f(O), 9 = -0.5(+) and 
9 = —l.O(n) respectively. The solid lines are obtained from the mean field solution Eq. ([T9| ) 
From Eq. ([□]) one finds: c curvatU rc(# = 0.5) = 0.67, c^ urvature (0 = 1-2) > 1, c^ urvature (0 = § ) = 
0.55, c curvature (0 = -1.0) = 0.23, and c curvature (0 = -0.5) = 0.27. 

large H the 3-dimensional nature of the liquid meniscus becomes more important and the 
2-dimensional mean field picture breaks down. 

In conclusion we have proposed a macroscopic Hamiltonian approach to the pinning of 
solid-liquid-vapour interfaces due to presence of stripes of obstacles. We find that curvature 
effects play a crucial role for the transition from the single pinning regime to the collective 
pinning regime. We have proposed an estimate of the density of pinning stripes for which the 
collective pinning happens, which is in good agreement with the simulations of the Langevin 
equation. We conjecture the same estimate to be valid in the case of random pinning sites. 
Finally a 2-dimensional mean field solution has been introduced which for small values of 
\9' — 9\,k gives excellent approximation for the height of the triple line. 

J.V.A. wish to acknowledge support from the Danish Natural Science Research Council 
under Grant No. 9400320. 
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